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^ ■ This paper reports on the detailed rp-process reaction flow on an accreting neutron star 
and the resulting ashes of a type I X-ray burst. It is obtained by coupling a 298 isotope 



reaction network to a self-consistent one-dimensional model calculation with a constant 
accretion rate of M = lO^^g/s. 
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d ■ 1. Introduction 

A neutron star in a low mass X-ray binary accretes a mixture of hydrogen and he- 
^ I hum from the secondary star. The photospheric impact ionizes the matter completely 
I afterwhich the matter spreads around the star, where it sinks down and develops a ther- 
monuclear instability after a few hours or days |1 . The instability causes an explosive 
runaway, where rapid and consequtive hydrogen and helium captures compete with j3~^- 
decays to fusion the accreted matter into heavier proton-rich material 

The specific reaction flow depends on the ignition composition, the temperature and 
density evolution and the capture-rates, decay-rates and photodisintegration-rates of the 
nuclei in the reaction flow [3 . However, near the dripline most of these rates are unknown 
which forces models to rely on theoretical predictions based on global models. 

In recent years new experimental developments - in particular the introduction of ra- 
dioactive beams - have made it possible to measure these rates experimenally |3]. Yet 
radioactive ion beam experiments take a long time to prepare and execute, so it is impor- 
tant to accurately establish the reaction flow to determine which reactions are critical in 
determining the reaction flow in the X-ray burst. 

The charged particle reactions of the rp-process are strongly temperature-dependent. 
Therefore realistic values of the temperature evolution during the burst must be provided. 
Currently, this is done best by a self-consistent one- dimensional model. 
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Figure 1. From left to right (solid line): 
y = 2.1 X lO^g/cm^ (surface), y = 9.5 x 
lOVcm^ y = 1.9 X Wg/cm^, y = 3.3 x 
lO^g/cm^, y = 6.2 x 10''g/cm^ (above igni- 
tion), y = 8.3 X lO'^g/cm^ (ignition point), 
and y = 1.1 x lO'^g/cm^ (ocean), where 
y = pdr is the column depth. The 
dashed line indicate the region which is 
convective during the rising of the burst. 



Figure 2. The mass fraction of the ashes 
about midway in between bursts plotted as 
a function of nucleon number: The peaks 
correspond to ^He; ^^C; ^^Si; ^^S and as- 
sociated isotopes reacting in this region; 
40Ca; 46Ti; and ^SNi, ^OZn, ^^Ge, ^^Se, ^^Kr, 
and ^^Sr and their respective daughters. 



2. The thermonuclear runaway, the reaction flow, and the resulting ashes 

This paper uses the same model as |5J. It considers a neutron star accreting at a 
constant rate of M = 1 x lO^^g/s which corresponds to the inferred accretion rate of 
GS 1826-24 jB]. Fig. Q shows the density and temperature for seven different depths as 
a function of time for a full cycle. ^ The mixed H/He bursters have the same kind of 
reaction flow throughout the atmosphere; where the maximum temperature at a given 
depth generally determines how much material is processed. It also determines the extent 
of the flow, except for deeper layers, where the exhaustion of hydrogen limits the flow to 
heavier isotopes. 

The composition is set by steady quiescent burning of freshly accreted hydrogen and 
residual hydrogen from the surface ashes from the previous burst steadily increasing 
the fraction of helium through the /3^-limited hot CNO cycle [7_ . At this temperature 
(T ~ 0.2GK c.f. Fig. ^ the ^^0(a, 7)^^Ne-rate dominates the otherwise much faster 
^^0(a,p)^''F-rate ^8J and it is responsible for generating a distribution of seed nuclei up 
to the well-bound ^°Ca isotope for the rp-process which are important to the runaway jU] 
that ensues at a column depth oiy = 8.3 x lO^g/cm^. 

^Densities and temperatures for different depths as a function of time along with initial compositions are 
available for post-processing calculations at ,http:/ /www.nucastrodata.org^ 
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Initially the breakout of the hot CNO cycle is led by ^^0(a,7) ^^Ne(p, 7) ^°Na(p, 7) 
^^Mg(p, 7)(7,p) ^^Al which is blocked by photodisintegration, whence the flow proceeds 
via ^^Mg(/3"'", Ti/2 = 0.122s) ^^Na(p, 7) ^^Mg. However, as the temperature rises and ^^O, 
which depends on the T1/2 = 76.4s ^^O beta-decay, is depleted, the reaction flow follows 
'^He(2Q;,7) ^^C(p, 7) -'^^N(p, 7) ^^0(q;,p) ^'^F{p,'~f) ^^Ne, which decays and returns via the 
hot CNO bicycle (see 0) until the temperature reaches T ~ 0.9GK at which point the 
flow joins the other breakout with ^^Ne(a,p) ^-^Na(p, 7) ^^Mg. 

From ^^Mg two processes lead the flow toward heavier nuclei, namely the (a,p)-process 
and the rp-process, whose different time-scales may lead to observable structure in the 
burst luminosity |5|lUj . The (a,p)-reactions primarily occur on longer-lived nuclei in the 
rp-process flow, but since the cross section decreases as the Coulomb barrier increases 
due to the increasing charge of the target as the flow proceeds, these reactions are only 
effective for A < AO. 

Higher temperatures make the proton-captures in the rp-process faster than the cor- 
responding /3^-decays and move the reaction flow towards the proton dripline. Yet for 
temperatures above ~ IGK the tail of the Boltzmann distributed photons becomes suf- 
ficently energetic to photodisintegrate the more weakly bound proton-rich nuclei, so the 
flow moves away from the dripline ^T]. If the flow branches into different paths, the 
branching ratio depends on the rate of the /3+-decay and the proton-capture, whence it 
becomes important to determine these rates accurately. 

Here the rp-process flow branches into either ^^Mg(p, 7) ^^Al(p, 7) ^^Si(/3^, T1/2 = 
0.102s) 24A1 or 22Mg(/?+,Ti/2 = 3.86s) 22Na(p,7) ^^Mg(p,7) ^^Al. This is followed by 
^■*Al(p, 7) ^^Si, which is in (p, 7)(7,p)-equilibrium with ^^P, where at high temperature, a 
fast additional proton capture creates the short-lived ^^S, which decays to ^''P. Otherwise 
the flow follows ^^Si(/3^, Ti/2 = 0.634s) ^^Al(p, 7) ^^Si, which may either decay or go to 
^^P after an additional proton capture. The flow continues with ^^Si(/5"'", T1/2 = 2.23s) 
26Al(p,7) ''Si(p,7) 28p or 27p(p,7) ^8S{(3+,T,/2 = 0.125s) ^sp or 27p(/5+,Ti/2 = 0.260s) 
^'''Si(p, 7) ^^P, which after a proton capture reaches ^^S(/5+,Ti/2 = 0.187s) ^^P(p, 7) ^°S. 

The decay half-time of ^^S is rather long (T1/2 = 1.18s), so initially the flow pro- 
ceeds from ^^Cl which is in (p, 7)(7,p)-equilibrium, since the Q- value is only 296 keV, 
until the temperature becomes sufficiently high to completely photodisintegrate most of 
the ^^Cl. At this point the flow must proceed through the sulfur waiting point i.e. 
(/3+,Ti/2 = 1.18s) 30p(p,^) 3is(p,7) 32ci(/?+,Ti/2 = 0.298s) 32S(p,7) 33ci(p,7) 34 Ar. 
The proton-capture Q- value on ^^Ar is even lower, {Q = 78 keV), so photodisintegration 
is stronger which makes ^"^Ar another waiting point. Therefore the flow proceeds via 
^^Ar(/5+, Ti/2 = 0.844s) ^''Cl(p, 7) ^^Ar(p, 7) ^^K. Here the flow branches and proceeds 
either via ^^K{f3+,Ti/2 = 0.342s) 36Ar(p,7) or ^6K(p,7) 37Ca(/3+, T1/2 = 0.175s) ^^K, 
which captures a proton to '^'^Ca. ^^Sc is almost proton unbound, so the flow must wait 
for 38Ca (Ti/2 = 0.440s) 3^K(p,7) ^^Ca. 

Here the flow branches to either ^^Ca(p, 7)(7,p) ''°Sc(p, 7) Ti (/?"'', T1/2 = 0.080s) 
4iSc(p,7) ^^Ti or 39Ca(/5+,Ti/2 = 0.860s) ^^K{p,-f) ^°Ca(p,7) ^^Sc(p,7) ^^Ti, where the 
latter path is active during the initial stage of the runaway burning the ^°Ca formed dur- 
ing the quiescent stage. Now the flow branches to either '^^Ti(/9"'", T1/2 = 0.199s) '^^Sc(p, 7) 
^3Ti(p, 7) ^^V or 42Ti(p, 7)(7,p) ^^V{p, a) ^^Cr(/3+, T1/2 = 0.053s) ^^V, where they connect 
and continue with ^^V(p,7) ^^Ci{p+,Ty2 = 0.050s) ^^V(p,7) ^^Cr, which decays to ^^Ti 
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via ^^V. 

The flow proceeds to ^^Mn via ^^Ti(p, 7) ^'^V(p, 7) ^^Cr(p, 7) "^^Mn or ^^Cr(p, 7) '^''Mn(p, 7) 
^8Fe(/?+, Ti/2 = 0.048s) ^^Mn, where it meets ^^V(p, 7) ^^Cr(p, 7) 48Mn(p, 7) ^^Fe{P+, Ty2 = 

0. 070.) ^^Mn. At this branching point the proton capture is weaker than the /3"'"-decay! 
Therefore the main flow proceeds via ^^Mn(/3"'"Ti/2 = 0.382s) ^^Cr(p, 7) ^°Mn(p, 7) ^^Fe{p, 7) 
52Co(p,7) s^Ni. 

At this point the flow becomes more dependent on the nuclear masses, which determine 
the proton-capture Q- values, since the photodisintegration of the compound nucleus deter- 
mines the specific waiting points. The waiting points are: ^^Ni, ^^Ni, and ^^Ni, where the 
flow either proceeds with two fast proton captures to ^"^Zn or awaits the T1/2 = 0.212s de- 
cay. Similiarly ^''Zn, ^^Zn, ^^Zn, and ^°Zn are waiting points. However, ^''Zn has a halflife 
of 2.38 minutes, so the flow must proceed with two proton captures reaching another row 
of waiting points: ^^Ge, ^^Ge, and ^^Ge which lives for 2.67 minutes. A proton-capture 
Q-value of 80keV essentially stops the reaction flow here, which most slowly decay out 
of ^"^Ge to form heavier nuclei. This means that most burst ashes are concentrated on 
the ^^Ge daughters as can be seen from fig. |21 Further waiting points which are reached 
during the ~ 150s cooling time of the burst are ^^Se, ^®Se, ^^Kr, and ^^Sr after which the 
flow becomes negligible. 

Ashes heavier than A = 40 are primarily made above the ignition point during a burst 
and subsequently accreted into the ocean, whereas lighter ashes are made directly in the 
ocean during the helium burning in the heat bath following a burst (see flg. |2I). 

3. Conclusion 

In this paper we have determined the rp-process reaction flow on an accreting neutron 
star using a self-consistent one-dimensional model. Waiting points and important reac- 
tions have been identifled. The ashes are concentrated around A ~ 30 and A = Q4 and 
show residual amounts of about 3% ''He and 3% ^^C. This disagrees with earlier calcula- 
tions of one-zone models which predicted much heavier ashes but lacked a self-consistent 
prescription for the heat transport. The present result corroborates other one-dimensinal 
models 
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